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The one body density matrix, momentum distribution, natural orbits and quasi hole states of O 
and 40 Ca are analyzed in the framework of the correlated basis function theory using state dependent 
correlations with central and tensor components. Fermi hypernetted chain integral equations and 
single operator chain approximation are employed to sum cluster diagrams at all orders. The optimal 
trial wave function is determined by means of the variational principle and the realistic Argonne Vg 
two-nucleon and Urbana IX three-nucleon interactions. The correlated le O momentum distribution 
is in good agreement with the variational Monte Carlo results and shows the well known enhancement 
at large momentum values with respect to the independent particle model. A similar behavior is 
found in 40 Ca. The relative importance of the different types of correlations (mainly Jastrow and 
tensor) on the momentum distribution appears to be similar in the nuclei and in nuclear matter. 
Diagonalization of the density matrix provides the natural orbits and their occupation numbers. 
Correlations deplete the occupation number of the first natural orbitals by more than 10%. The 
first following ones result instead occupied by a few percent. The single particle overlap functions 
and the spectroscopic factors are computed in the correlated model for both nuclei and compared 
with previous estimates. Jastrow correlations lower the spectroscopic factors of the valence states 
by a few percent (~ 1 — 3%) with respect to unity. An additional ~ 8 — 12% depletion is provided 
by spin-isospin tensor correlations. It is confirmed that a variational teatment of short range 
correlations does not explain the spectroscopic factors extracted from (e, e'p) experiments. Such 
an approach corresponds to the zeroth order of the correlated basis function theory and two-hole 
one-particle perturbative corrections in the correlated basis are expected to provide most of the 
remaining strength, as in nuclear matter. 



I. INTRODUCTION 



The notion of nuclei as a set of mutually, strongly interacting particles is now generally accepted and it is widely 
recognized that correlations beyond the mean field play a substantial, if not decisive, role in the microscopic description 
of nuclear properties. This is intuitive when energies are evaluated, but there are also clear signatures of the presence 
of correlations in some quantities related to the behavior of the single nucleon in the medium. For instance, it is 
known that the one nucleon momentum distribution (MD) has dominant high momentum components which are due 
to short range nucleon nucleon (NN) correlations and are not describable in any independent particle model (IPM) 
[0,1). In addition, NN correlations may be responsible for the reduction of the spectroscopic strengths of the hole 
states ||. 

The ideal way of describing the nucleus would consist in solving the many-body Schrodinger equation with realistic 
interactions. Exact solutions have been obtained in light nuclei up to A=8 within a variety of methods: quantum 
Monte Carlo Faddeev-like |(| and correlated hyperspherical harmonics expansions. 

The treatment of heavier nuclei has not yet attained the same degree of accuracy as the light ones. However, the 
situation is rapidly improving, at least for doubly closed shell nuclei. The 16 nucleus has been studied by variational 
Monte Carlo || (VMC) and coupled cluster Q methods. The variational theory underlying the correlated basis 
function (CB F) jlO(] method has been used in doubly closed shell nuclei by Fermi hypernetted chain (FHNC) integral 
equations JTT| , [l2| . The same accuracy as in the best variational studies of nuclear matter has been obtained in 16 O 
and 40 Ca by using spin and isospin dependent correlations and modern, microscopic potentials flS] ]. To this aim, the 
nuclear matter single operator chain approximation J14| has been extended to finite nuclear systems to deal with state 
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dependent correlations in the framework of the FHNC technique M] (FHNC/SOC). The variational and the coupled 
cluster approaches have provided an overall satisfying description of such ground state properties as binding energies, 
one- and two-body densities and structure functions. 

In medium heavy nuclei the microscopic approach has to be compared with the independent particle model and its 
description of various nuclear features. The use of this model is so wide and well-established that great part of the 
nomenclature on medium and heavy nuclei is based upon concepts defined within the IPM itself. For example, the 
ideas of collective excitations, single particle levels, occupation probabilities and spectroscopic factors are meaningful 
only in a theoretical framework where the IPM is a first order approximation. 

The sources of effects beyond the IPM description are generally termed as correlations. This name is hiding two 
rather different kinds of physical effects. Long range correlations are generated by the so-called residual interaction, 
neglected in the IPM and reponsible for the collective excitations. The short range correlations (SRC) are mainly 
produced by the repulsive core of the nucleon-nucleon interaction. Moreover, important non central components at 
intermediate range arc needed for most realistic descriptions. While collective phenomena have been known since the 
early times of nuclear physics, effects produced by the SRC are more difficult to be experimentally singled out. Only 
in these last few years, thanks to the advances of the experimental techniques, there have been some consistent effort 
aimed to their identification. 

The sensitivity of the two-nucleon emission cross section to the SRC is evident Jll|, but also other quantities 
related to the behavior of the single nucleon in the medium seem to depend on them. For instance, (e, e'p) data 
in the quasi-elastic region need a consistent reduction of the IPM hole strength to be reproduced JL7]|. The same 
holds for the electromagnetic form factors of low lying states with high angular momentum pH| |. Furthermore, charge 
density distributions obtained by elastic electron scattering experiments are, in the nuclear interior, smaller than 
those predicted by the IPM [|l9| . All these facts could be explained by assuming occupation probabilities of the single 
particle levels different from that of the IPM J2(J . 

From the theoretical point of view, the basic quantity to be investigated in order to verify the hypothesis of partial 
occupation probability is the one body density matrix (OBDM), p(ri,ri/). defined as: 

p(n, r r ) = (tfoWlatfnMriOltfotA)) , (I) 

where ^o(A) is the ground state A-body wave function and a^(ri) is the creation operator of a nucleon at the 
position ri . Much theoretical effort has been devoted to understand the behavior of the OBDM and of the momentum 
distribution (MD), n(k), given by the Fourier transform of p(ri,ry). The MD is also obtained by the energy integral 
of the spectral function, that is often used in plane wave impulse approximation to study inclusive and exclusive 
reactions. In some approximations and kinematics, the MD is directly employed. The natural orbits (NO) plj , with 
their occupation numbers (n Q ), are defined as the basis where the OBDM is diagonal. In the IPM, the nuclear ground 
state is described by a Slater determinant of fully occupied single particle (SP) wave functions below the Fermi surface, 
ap. In this case, the NO and the SP w.f. coincide and n a < aF = \ and n Q>QF =0. Deviations from this situation are a 
measure of the correlations, since they allow higher NO to become populated with n a ^0. 

Quantities not directly related to the OBDM, but accessible to the experimental investigation by means of knock 
out experiments, are the quasihole (QH) wave functions, iph(r), defined as the overlaps between and the hole 
states, i&h, produced by removing a nucleon from the position r. From (e, e'p) experiments it is possible to obtain 
an accurate determination of the QH overlap functions |p2| . Their normalizations give the spectroscopic factors, Sh, 
that deviate from unity (the IPM value) because of various effects, from center of mass to correlation corrections. 
Typical values extracted from the experiments are Sh ~ 0.6 — 0.7 |p3| , and their correct calculation represents a severe 
challenge for the available many-body theories. 

The FHNC theory for the OBDM in correlated infinite matter was developed by Fantoni |24f. Subsequentely, 



25| ] and the spectral function |26J . The technique 
1 nuclei described by a correlated wave function 



FHNC/SOC was applied to nuclear matter to calculate both the MD 
has been then extended to evaluate the OBDM in doubly closed she ' 
containing scalar, spin-isospin independent (or Jastrow) correlations [{27|| . Isospin dependence was later introduced 
to treat closed shell nuclei in the jj coupling scheme |IJ]. The 16 O momentum distribution has been calculated by 
using state dependent correlations within the VMC approach. Low order cluster expansions have been employed 



in closed shell nuclei 51 29j] and in sp and sd shell nuclei with N=Z []30|. Local density approximation was used to 
estimate the MD in medium-heavy nuclei J3l| from the nuclear matter results. The effect of short range correlations 
on the OBDM and the NO in 16 O has been studied by means of the Green function method p2[ . As far as the overlap 
functions are concerned, in the Li(e, e'p) 6 He reaction the experimental data have been succesfully compared with the 
VMC, parameter-free theoretical predictions 1 33 1 . VMC has been used also in the analysis of the 16 0(e, e'p) knockout 



experiments |34| . A relationship between the OBDM and the QH overlaps |35| has been recently exploited [p6|]37| l to 
evaluate the spectroscopic factors for the same reaction, using several correlated models for the density matrix. The 
effect of long-range correlations has been considered in Refs. p8,pq]. 
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In this paper we extend the FHNC/SOC methodology to deal with the OBDM of 16 and 40 Ca using wave functions 
with central and tensor correlations. The correlated A-body wave function, ^(l^ 2. ..A), is written as: 



tf (l,2...A) = G(l,2...A)$o(l,2... J 4). 



(2) 



where G(1,2...A) is a many-body correlation operator acting on the mean field wave function $o(l, 2. ..A), given by 
a Slater determinant of single particle wave functions, <fi a {i)- The correlation operator is given by a symmetrized 
product of two-body correlation operators, Fy, 



G(1,2...A) =S 



In the most sophisticated variational calculations (both in nuclear matter and in nuclei) Fij assumes the form: 



F, 



E 



nnj)oi. 



(3) 



(4) 



where 



and Si 



O 



[l,o- l ■ (Tj,Sij, (L ■ S) 



(5) 



(3?, 



(Tj — (Ti ■ (Tj) is the tensor operator. The spin-orbit components of the correlation (p =7, 



have been omitted in this paper (f$ model), as well as spin-orbit and Coulomb interaction terms in the mean field 
potential generating the single particle wave functions, <p a {i). 

The variational principle allows for fixing the correlation functions, f p (r), and the single particle wave functions by 
minimizing the ground state energy, 



(*ol#l*o) 



(*o|*o) 

In our calculation we adopt a non relativistic nuclear hamiltonian of the form: 



H 



~2m ^r 1 



v » 



i<j<k 



(6) 



(J) 



Modern two-nucleon potentials, uy, are built on the very high quality phase-shift analyses of the NN scattering 
data 39 40[|. We have used the v' 8 reduction of the Argonne vi% |41| potential. For the three-nucleon interaction the 
Urbana IX model ||] has been adopted. 

The plan of the paper is as follows: section 2 is devoted to a brief description of the FHNC/SOC theory for 
the OBDM; in section 3 the results obtained for the 16 O and 40 Ca momentum distributions and natural orbits are 
presented and discussed; in section 4 we describe the FHNC theory for the overlap functions and give the results for 
the spectroscopic factors; the conclusions are drawn in section 5. 



II. FHNC/SOC THEORY FOR THE ONE BODY DENSITY MATRIX 

The one body density matrix ([!]) may be written as 

p(r 1 ,r 1 ,) = ^ J d 3 r 2 ...J d\ A ^ (1, 2, ..v4)* (l', 2, ..A) , 



where JV = J d 3 r\... j d 3 rvi| , I'o| 2 - The one body density, pi(ri), is the diagonal part of the OBDM, whose Fourier 
transform gives the momentum distribution, 



n(k) = 1 J d 3 n J d 3 r v p(r 1 ,ri,)e ik -( ri - r ^ . 
In the independent particle model, the OBDM is given by: 



(9) 
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pjpjtffa.n,) = 4(1)^(1') = (E^r(l)x-(l'))^o(ri,r 1 ,) ) (10) 

a err 

and 

Pi,/PA/(ri) - £ I^WI 2 = ^o(ri, ri ) , (11) 

a 

where Xo-r(l) is the spin-isospin single particle eigenfunction and v the degeneracy number (y — 4, for the doubly 
closed shell, N=Z nuclei we are considering). The second equalities in the above equations are valid since we are 
working in Is coupling. 

In analogy, we define the function iV(ri, ry) via the relation: 

p(r u r v ) = (^ fX lr(X)XarO'))N(T l ,T V ). (12) 
err 

We have already presented in Ref. |27| the FHNC theory of the OBDM for Jastrow correlated wave functions, 

Gj{l,2...A)=J[f 1 {r ij ), (13) 

having only the p = 1 scalar component in the correlation factor Fij. The density matrix is expanded in powers 
of the dynamical correlations, h(r) — [/ 1 (r)] 2 — 1 (/i-bond) and ui(r) — / 1 (r) — 1 (w-bond), and of the statistical 
correlations, N (ri,rj) (exchange or e-bonds). The expansion generates cluster terms classified according to the 
number of particles and to the number of the correlations. The FHNC equations allow for summing cluster terms at 
all orders. Details of the finite systems FHNC theory may be found in |27j and in Ref. JllJ. 

The FHNC/SOC equations for the more general f$ correlation were derived in Ref. ]15[ for the one and two body 
densities. In the case of the OBDM, we can write 

N(r u r r ) = £(n)[l + AC(ri)]£(ry)[l + A^iv)^' 1 '* 1 ' 5 (14) 
x [Noir^ry) - NZ eUc (ri,*i>)] + C(^)C(rr)e N ^' r ^ 

p>2 

N% u and iV£ c(ilc are nodal functions (see [^J) of the luuj and uj c lu c type, Q and A£° p are the central vertex corrections 
and their operatorial contributions [£5) and the matrices A? and A p are defined in (l^] . The components will be often 
labelled as c (p = 1) and a (spin), r (isospin) and t (tensor). 
The nodal funtions are obtained by the equations: 

tt.1')=EE / ^ a ^X2 s (r 1 ,r 2 )^(r 2 )[^ /(i ,(r 2) r 1 ,) + ^ tl> (r 2 ,r 1 0] , (15) 



xx' qr 



and 



qr J 
qr>1 J 

+ N^vyvy) + N'P^vr) + N^Jr^ry) . 

The index x(x') may assume the values x(x') — d,e and denotes the exchange pattern and the type of correlation at a 
specific point. In general, we may have: d— and w-vertices, if the point is reached by h— and w-bonds, respectively, and 
does not belong to any exchange loop; an e-vertex if the point belongs to a closed exchange loop and it is reached by 
two e-bonds; c- and w c -vertices if the point belongs to an open exchange loop and it is reached by a single exchange 
line. The allowed (xx ) combinations are: dd, de, ed. We indicate with V^, the vertex corrections, and with the 
angular couplings. The expressions of these equations are all given in [19]. 
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The above equations are derived in the FHNC/0 approximation, which does not include the contribution of the 
elementary diagrams. A detailed discussion of the importance of these diagrams and of the accuracy of the approx- 
imations used in solving the FHNC equations can be found in Refs. p^ , p3[ . The FHNC/0 choice has been used in 
studies of the equation of state of nuclear matter [fiof , where it was found that the elementary diagrams contribution 
is not important because of the relatively low densities of the system. In finite nuclei elementary dia gram s may play 
some role in the evaluation of the expectation value of potentials having strong exchange components . In general, 
a measure of the relevance of the missing diagrams is provided by the accuracy of the sum rules of the one and two 
body densities. 

The partial nodal functions, N^ x , are solutions of the integral equations: 



«ri,r 2 ) = EE / d ^&ly{ri^W^) [^ x (r 3 ,r 2 ) + A^ x (r 3 , r 2 )] 

or J 



(17) 



yy> qr 



The X£f 2 links are 



Xlf (n,r 2 ) = /£(!•!, r 2 )/£(n,r 2 ) - J\ft,(ri,r 3 ) , (18) 

X*AT 1 ,T i )=hl(T U T 2 )[}%(T U T 2 )NZ^ (19) 

/£(n,r 2 ) = P(r 12 ) + / c (r 12 )A^(n,r 2 ) , (20) 

and /£(n,r 2 ) = expjV£ d (ri, r 2 )]. Moreover, the N(X)P y (r 1 ,r 2 ) = N(X)P x (r 2 , r x ) property holds. The central X° x 
links are defined in p7[ , and the X? y and N£ y functions in |15| . 

For the w c c-type nodals we have N^^ = NP, cX + A/£ cp + Nf yx + N* p . and iV£ cP are solutions of 

NP cX (r 1: r 2 ) = J2 [ dV3^1c(ri,r 3 )^ r (r3)[^ c c c (r3,r 2 ) + iV^(r3,r 2 ) + 7V; :c (r3,r 2 )] A r (21) 

qr J 

+ E / ^i^A^r!,^^^^^ , 

qr>1 J 

A^ p (r 1; r 2 ) = W d'rs^X^Jv^V^irs) [-N (r 3 , r 2 ) + N* p (r 3 , r 2 ) + N c pp (r 3 , r 2 )] A r (22) 

qr J 

+ E / d 3 ^A^ cC (r 1 ,r 3 )K 9 c r (r3)[^(r 3 ,r 2 )+iV; p (r3,r 2 )] . 

The equations for N^ x and are given in |l5{ ]. The links X"^ are 

X^ c 2 (r 1) r 2 ) = / l S(r 1 ,r 2 )^(r 1 ,r 2 ) [NZJt u r 2 ) - N^, r 2 )] (23) 
+ [/ c (r 12 )^(r 1 ,r 2 )-l]^ cC (r 1 ,r 2 ). 

Again, X^ c is defined in fl27||. 



The vertex corrections, are discussed in 27 and A£° p is given by: 



AC P (ri) = C(ri) + Ey / ^r 2 {Xl d {v u v 2 )Nl d {v u v 2 )pl{v 2 ) (24) 



p>2 

+ K d ( ri , r 2 )C(n, r 2 ) + XP e ( ri , r 2 )^( ri , r 2 )] C d (r 2 )} 



/?f (the Jastrow part of the one body density) and Cd are given in 15|, and U° p is obtained by eq.(2.12) of [|27j with 

■ d 



the substitutions p\ — ► pi — Pi and —> CdU^. 
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III. RESULTS FOR THE MOMENTUM DISTRIBUTION AND NATURAL ORBITS 



In our work we used the Argonne v' s NN potential. This model is based upon the Dig potential and it is constructed 
by considering only the first eight operator terms, up to the spin-orbit ones. It reproduces the isoscalar part of the 
full interaction, v\g, in the S, P and 3 D\ waves and the 3 Di~ 3 Si coupling. Argonne v' 8 was introduced in Ref. || 
because its parametrization (simpler than that of other realistic potentials, since L 2 and (L ■ S) 2 components are 
missing) allowed for a large simplification in the numerically involved quantum Monte Carlo calculations. The v' s 
potential is slightly more attractive than wis, and, for this reason, the strength of the associated repulsive part of the 
three-nucleon force (Urbana IX model) has been increased by 30% with respect to the original version. The results 
presented in this paper have been obtained with this hamiltonian (A8'+UIX model). 

The correlation functions, f p {r), and the single particle functions, 4> a {i), are the two ingredients necessary to 
construct the many-body wave function (|^). We use a /g correlation, therefore, with respect to the structure of 
the hamiltonian, we neglect the spin-orbit components. The correlation is determined by minimizing the nuclear 
matter energy at the lowest order of the cluster expansion, considering the Fermi momentum, Uf, as one of the 
variational parameters. The resulting two-body Euler equations are solved with the healing conditions / (r > d±) = 1, 
f p>1 (r > dp) = 0, and requiring that the first derivatives vanish at r — d p . Only two healing distances are introduced, 
d c for the four central channels and c? t for the tensor ones, and they are variationally fixed. More details on this 
procedure are given in Ref. fl4f| for nuclear matter and in [ fl5[ for nuclei. 

The single particle wave functions have been obtained by solving the single particle Schrodinger equation with a 
Woods-Saxon potential, 

Vws(r) = — J° 5 w 1 • (25) 

l + exp[(r-i? )/ao] 

A full minimization for the A8'+UIX model has been obtained in Ref. |D| and it has provided a binding energy per 
nucleon, B/A, of 5.48 MeV in 16 and 6.97 MeV in 40 Ca, to be compared with the experimental values of 7.97 MeV 
( 16 0) and 8.55 MeV ( 40 Ca). These differences are comparable with those obtained in nuclear matter at the empirical 
saturation density, pnm =0.16 fm~ 3 , with the same hamiltonian. In fact, the FHNC/SOC nuclear matter energy per 
nucleon, E^m, is Bjvm =-10.9 MeV p3[ , against the empirical value of -16 MeV. 

The nuclear root mean square radii (rms) were found to be 2.83 fm in 16 and 3.66 fm in 40 Ca (the experimental 
radii are 2.73 fm and 3.48 fm, respectively). However, the one body densities at the variational minimum did not 
show a satisfactory agreement with the experimental ones. Moreover, shallow minima with respect to variations of the 
mean field parameters around the minimum itself were found in Ref. In particular, if one chooses a set of single 
particle wave functions which reproduces at best the empirical densities, the A8'+UIX model provided B/A=5.41 
MeV in 16 and B/A=6.64 MeV in 40 Ca, with rms( 16 0)=2.67 fm and rms( 40 Ca)=3.39 fm. Therefore, the density 
description has largely improved while the energy variations are kept within the accuracy of the FHNC/SOC scheme. 
The results presented in this paper have been obtained by means of this type of wave function, whose parameters are 
given in Table V of Ref. @. 

The one body densities generated by the FHNC/SOC scheme are shown in Fig. [|, where the solid lines give the 
densities obtained with the full correlation, the dot-dashed lines are those obtained with the Jastrow correlation (re- 
taining only the p =1 component) and the dashed lines are the IPM densities. The effect of the operatorial correlation 
is large with respect to the Jastrow case, that is hardly distinguishable from the IPM one. The comparison with the 
experimental results has been presented in Ref. O], where the proton densities are folded with the electromagnetic 
nucleon form factor. 

The momentum distributions are given in Fig. ^. Again the solid and dot-dashed lines are the fully correlated and 
Jastrow results, respectively, while the IPM ones are shown as dashed lines. The squares are the VMC results || for 
16 O obtained with the Argonne vn |Q| NN interaction. 

The MD is normalized as: 



1 = j^rJ d kn{k) > (26) 

where v is the spin-isospin degeneracy. For the Jastrow correlation this normalization is satisfied within the 0.2-0.3 
% in both nuclei, while for the fe model the error is ~3 % in 16 and ~2 % in 40 Ca, reflecting the approximations 
of the SOC approach. 

In Fig. |^ the MD of 16 O (thin continuous line) and that of 40 Ca (thin dashed line) are compared with those of 
nuclear matter normalized as in ( p6|) and calculated in FHNC/SOC framework by using the same interaction. It is 
worth noticing that the differences between the Jastrow and the fe correlations are similar in the infinite and finite 
systems and that the three cases show an analogous behavior at large momentum values. This momentum region is 
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dominated by the short range structure of the nuclear wave function, which is heavily affected by the NN correlations. 
The effect appears to be, to a large extent, independent on the nucleus. A similar behavior was found in Ref. [0, 
where the comparison was made among the 4 He, 16 O and nuclear matter momentum distributions. With respect to 
the Jastrow estimates the non central, tensor correlations enhance the tails of the MDs by a factor 2-3 slightly smaller 
than the one found in Ref. [||, which is roughly ~4. The difference may be understood in terms of the stronger tensor 
force of the Argonne vu potential adopted in that reference. Part of the discrepancy may also be ascribed to the 
presence of spin-orbit correlations in the wave function of Ref. However, we notice that, in the same paper, it 
was found that these correlations contribute to the kinetic energy by only ~ 1% (0.4 MeV/A out of a total kinetic 
energy of 34.4 MeV/A). 

A more demanding sum rule for n(k) than that expressed by eq.(p6|) can be obtained from the kinetic energy, T. 
The kinetic energy per particle can be evaluated via the MD as: 



t h 2 



A 2m (2tt) 



d 3 kk 2 n(k) = T MD - (27) 



The value of T/A can also be computed in the FHNC/SOC framework (Tfhnc), for example by means of the Jackson- 
Feenberg identity, as it has been done in Ref. for the A8'+UIX model. The differences between Tmd and Tfhnc 
are a severe measure of the importance of the approximations made in the cluster expansion (FHNC/0 and SOC). 
For the Jastrow cases the relative disagreement, ST = \Tmd — Tfhnc\/TfhnCi is < 5% ( Tmd=20.52 MeV and 
T FH jvc=19.57 MeV in 16 and T MD =22M MeV and T F hnc=22.05 MeV in 40 Ca) and it is due to the absence of the 
elementary diagrams in the FHNC/0 truncation. In the f 6 model we obtained Xm,d=29.42 MeV, Tfhnc— 32.64 MeV 
and ST =9% in 16 and T MjD =36.63 MeV, T FHNC =38.15 MeV and ST =9.6% in 40 Ca. This larger disagreement is 
due to the SOC approximation. The Tmd value is largely influenced by the behavior of the momentum distribution 
at high fc-values. The contribution of the k >5 fm _1 tail has been evaluated by an exponential extrapolation of the 
computed MD. The tail contributions for the Jastrow and operatorial correlations are about 5% and 10% of the total 
Tmd, respectively. So, we believe that the uncertainty in Tmd related to the MD tails may be fixed to a few percent 
in both cases. As a additional check of the numerical accuracy of the algorithm used to evaluate the momentum 
distribution, we have verified that Tmd coincides with Tfhnc for the IPM. 

The NO and their occupation numbers are obtained by diagonalizing the OBDM: 

Pl (r x , r y) = J2 n*<£° (r i ) f cp N a ° (r y) , (28) 

a 

We treat spherical nuclei in Is single particle coupling, saturated in both spin and isospin. For this reason the 
spin-isospin part of eq. (|i~2| ) provides the degeneracy i/=4. Because of the spherical symmetry the function N(r\, ry) 
of eq. ( |l2| ) can be expanded in multipoles, and we obtain for the OBDM: 

21 + f 

Pl {vyv v ) = ^——i^cosfluOpKn.ri/) , (29) 
i 

where Pi(x) represents the Legendre polinomials and Q\y is the angle between ri and vy. 
Exploiting again the spherical symmetry, the natural orbitals can be written as: 

^2 nlm (r) = €i°(r)Yi m (r)x,r, (30) 

where we indicate with Y; m (r) the spherical harmonics and with \ aT the spin-isospin part of the wave function. The 
normalization condition is: 

1= \r 2 dr\^°{r)\ 2 . (31) 



Therefore we obtain: 



Pl {ryry) = vJ2 n m €l° (n) • (3 2 ) 



The nZ-natural orbitals and their occupations have been obtained by discretizing and diagonalizing the matrix 
Pi{ r i7 r v) m a 100 x 100 equally spaced grid, up to r max = 6(7) fm for 16 ( 40 Ca). 

The first three NO of 16 and 40 Ca, calculated for the three lowest Z-values are shown in Figs. [I] and [|. In Table || 
the occupation numbers of the various NO for Jastrow and / 6 correlations are presented. 
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The effect of the correlations on the shell model orbitals are mainly visible in the short range part of the Is state, 
making this NO more localized than its IPM counterpart. The shape of the other IPM states is barely influenced 
by the correlations. The occupation of the NO corresponding to the shell model ones is depleted by as much as 22% 
(the 2s state in 40 Ca). In contrast, the mean field unoccupied states become sizeably populated. The two effects are 
largely due to the tensor part of the correlation operator. 

It was pointed out in Ref. Jl5|] , that when there is more than one occupied state in the IPM for a given Z -value (as 
for the s-states in 40 Ca), then the natural orbitals may be qualitatively different from the IPM ones. In fact, any 
orthogonal combination of the mean field orbitals does not change the corresponding density matrix. New mean field 
s-orbitals, 4> ns , consistent with pipm(ti, r i')> can be obtained in 40 Ca by the transformations: 

4>u (r) = cos(a)^i s (r) + sin(a)</> 2s (r) , (33) 



4>2s (r) = cos(a)0 2s (r) - sin(a)0i s (r) , (34) 

for any choice of the angle a. The IPM orbitals of the 1=0 panel in Fig. || are obtained by a numerical diagonalization 
of the one-body density matrix and roughly correspond to cos(a) = .8. 

The 16 O natural orbits have been evaluated in Ref. |3^] within a Green function approach and using the one- 
boson-exchange Bonn B potential of Ref. jf6|. The authors find the n =1 NO more populated than the CBF ones for 
the occupied states in the shell model approach (I =0, 1), and, consequently, lower occupations for all the remaining 
orbitals. Specifically, the n=l Green function results are: ni s =.921, n lp =.941 and The lp (Id) occupation 

number has been taken as the average of the lpi/2 and IP3/2 (1^3/2 and 1^5/2) orbitals given in the reference. The 
discrepancies are probably to be ascribed more to the different potentials adopted, rather than to the methodologies. 
The A8'+UIX model induces stronger correlation, so giving a larger depletion of the lowest NO. This effect was also 
found in the study of 3 He atomic drops of Ref. f45|| , where the strong repulsive interaction between the 3 He atoms 
depletes the shell model occupations by 15-46 %. The CBF total occupation numbers in the Z— tli orbitals for 16 and 
40 Ca with different correlations (fe, f^, without tensor components, and Jastrow, J) are given in Table [n| together 
with the 16 O Green function ones from Ref. j32|. It appears clearly that the longer ranged tensor correlations are 
responsible of most of the deoccupation of the shell model natural orbitals in favour of the higher ones. 

Fig. ^ presents a comparison between the 16 NO in FHNC/SOC and in the lowest order approximation. This 
approach consists in truncating the cluster expansion at the first order in the dynamical correlation lines p7| , p8| , p9[ 
and it has achieved a certain degree of popularity because of its simplicity. The approximation provides a good 
description of the n=l NO corresponding to the occupied shell model states, but it fails to reproduce the other ones. 

To conclude this section, we give in Fig. fj] the partial wave decomposition of the correlated and IPM one body 
densities in terms of their natural orbits. 



IV. QUASI HOLE STATES AND SINGLE PARTICLE OVERLAP FUNCTIONS 

A considerable amount of information on the properties of the single nucleon in the nuclear medium can be deduced 
from (e,e'p) reactions. These experiments have been analyzed to extract the quasiholc function, iph (r), given by the 
overlap between the A-body ground-state and the (A-l)-body hole state of the residual system. 

In a fixed center reference frame, the QH function is defined as: 

In doubly closed shell nuclei in the Is coupling scheme it is possible to separate the radial dependence of the QH 
function from the angular, spin and isospin ones, as: 

ip h {x) = i/Jh{r)Y lm (f)xcTT = i>h(r)yim<7r(r) ■ (36) 

In the IPM, the QH overlaps are simply the shell model functions and ip I h PM (r) = Rh=ni{r), where Rh(i~i) is the 
radial part of <fi a (i). 

In CBF theory $ (1,2...A) is given by (g) and 

2...A - 1) = G(l, 2...A - l)$ h (l, 2... A - 1), (37) 

where 2.. .A — 1) is a Slater determinant obtained by removing from $o(l, 2. ..A) a nucleon in the state h. 
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In order to develop a cluster expansion for ^ph{r) it is convenient to rearrange ( p5| ) as: 

Mr) = X h (r)Ml' 2 , (38) 



where 



and 



^ = VI PlNni ' (39) 



^ - ' (40) 

Cluster expansions are used to compute Xh and A4, along the lines followed in Ref. |5(J to evaluate the overlap 
matrix elements in the CBF approach to the nuclear matter spectral function. 

The expansion for X% is linked, in the sense that disconnected diagrams in the numerator are exactly canceled by 
those coming from the denominator. Its FHNC/0 expression, when only Jastrow correlations are considered, is: 

*»i(»0 = CV) Um(r) + J cPnRMPiicoae) [<7« d (r, ri)C^(ri) (41) 

where 8 is the angle betwen r and ri, and p IPM is the A-l one-body density matrix in the independent particle 
model, 

^ PM (ri,riO = E^( 1 )^( 1 ')- ( 42 ) 

In J\fh only those diagrams from the denominator containing explicitely the /i-orbital survive. The Jastrow, FHNC/0 
expression is: 

[Aft]" 1 - J d?rCf{v) ||^(r)| 2 + J dVi^r^Cn) [gf d (r, T^Cf^) (43) 
x (-pf PM (r, n) + N^(r, n)) + pf PM (r, n) - A^(r, n) - JV#(r, n)] } . 

The FHNC quantities entering X^ and A"^ correspond to those given in Refs. |2^JTl[], but evaluated with the A-l 
densities, p ipm {y u y v ) and Pi. IPM ( r i) = J2 a ^h I'M 1 )! 2 - In absence of correlations, X^{r) -> R n i{r) and A^ -> 1. 

The FHNC/SOC expressions of A"„; and A^z, with operatorial correlations, are given in the Appendix. 

The quasihole normalization gives the spectroscopic factor, Sh, 

S h = J r 2 dr^l{r). (44) 

In a fixed center IPM (as the one we adopt as model function), S^ PM = 1. Center of mass (cm) corrections are 
sources of deviation. In the harmonic oscillator model they enhance Sh for the valence hole states (those with the 
largest oscillator quantum number, N v ) by a [A/(A — 1)]^" factor [^9|. As a consequence, the cm-corrected Ip- 
shell spectroscopic factor of 16 is S^ p ° cm — 16/15 ~ 1.07, while the average between the 2s and Id states in 40 Ca 
^ s ^2sfid cm = (40/39) 2 r-j 1.05. More realistic Woods-Saxon orbitals do not allow for an analytical treatment of 
cm effects, which have to be computed numerically. It has been found that in 16 the lp WS spectroscopic factor 
practically coincides with the HO one [B4[ 



The correlated spectroscopic factors (without cm corrections) in 16 O and 40 Ca are given in Table III for /g, fi and 
Jastrow correlation factors. Jastrow correlations marginally reduce Sh (at most 3%). The Jastrow l-th spectroscopic 
factors may result slightly larger than the total occupation of the corresponding natural orbits, given in Table ||. 
A similar feature was found in Ref. |4g] . The small deviations of the Jastrow model from the correct behavior are, 
in our opinion, well within the accuracy of the numerical procedures we have adopted and of the approximations 
in the cluster summation. Central spin-isospin correlations (/4 model) also provide a few percent depletion in the 
valence states, whereas the tensor ones (fe) give most of the reduction bringing Si p in 16 to 0.90 and S2 S and Sid 







in 40 Ca to 0.86 and 0.87, respectively. The lp CBF 16 result is in complete agreement with the VMC estimate [ [34[ . 
The influence of the operatorial correlations is much larger in the low lying states, whose spectroscopic factors are 
drastically reduced by both central and tensor components: Si s in 16 is 0.70, S± p and Si s in 40 Ca are 0.58 and 
0.55, respectively. An analogous behavior was found by Benhar |j5l| , who first used low order cluster expansions to 
estimate Sh in the 12 C(e, e'p) 11 B reaction with state dependent correlations and found Si p =0.55 and 518=0.25. 

Results similar to those presented in this paper have been obtained in Refs. J3^^], where the 16 Si p has been 
extracted by several models of OBDM p5| . In particular, in both references it is confirmed that correlation effects 
on this spectroscopic factor are dominated by tensor components. The lowest order truncation of the OBDM cluster 
expansion adopted in Ref. |37]] provides Si Pi lo ~ 0.90, in agreement with the FHNC/SOC results. However, the 
authors find S\ S: lo ~ 0.86, in contrast with the 0.70 FHNC/SOC value. The origin of this large difference may lay 
in the lowest order approximation in the cluster expansion. This issue is presently object of investigation. 

The latest experimental extraction of S p from the 16 0(e, e'p) 15 N reaction |52| reports S'p 1 . 2 =0.61 for the l/2~ 
ground state in 15 N and S Pa/2 (6.32)=0.53 for the lowest 3/2~ state at 6.32 MeV. This state exhibits 87% of the total 
S P3/2 strength, that is fragmented over three states at 6.32, 9.93 and 10.70 MeV. So, the total S P3/2 may be estimated 
to be S P3/2 =0.53/0.87=0.61 @. 

A corresponding situation is met in the 40 Ca(e, e'p) 39 K reaction p3|] , where the transition to the ld^M ground state 
gives Sd 3/2 ~0.61±0.07, while the FHNC/SOC value is Sd =0.87. The 40 Ca spectroscopic factors have been computed 
by the low order cluster expansion of the OBDM in a Jastrow model in Ref. [p3| . The results are consistenly lower 



than the FHNC ones, reported in the 5th column of Table III. For instance, S2s,lo^0.95 and Si ( i.lo~0.91, whereas 



S2s,fhnc=0-98 and Sid,FHNC = 0.97 . The discrepancies are probably to be ascribed to the approximation used in 
the reference to evaluate the OBDM. 

The squared quasihole functions are shown in Fig.^. The solid and dot-dashed lines give the full fe and Jastrow 
results, respectively. The IPM estimates are given as dashed lines. The spin-isospin dependent correlations are the 
main responsible for the quenching of the IPM QH functions and, consequently, of the spectroscopic factor. The 
Jastrow components have little effect on the overlaps, and mostly in the valence states. In Ref. Jsij a Woods-Saxon 
potential was used to generate a single-particle wave function to fit the 16 0(e, e'p) 15 N cross section to the 6.32 
MeV state with (6.32)=0.53. IV'ws'l 2 is shown in the |^i p | 2 — 16 panel as stars. In order to give a meaningful 
comparison, we rescale \4>i p ,fhnc\ 2 by the factor 0.53/0.90. The result is shown as a dot-dashed line and it is in nice 
agreement with the empirical estimate. 

The knowledge of iph{r) may give access to the cross sections. However, both Coulomb distortion and final state 
interactions should be properly accounted for, by evaluating the Fourier transorm of a distorted overlap |j55| , to 
perform a quantitative comparison with the experiments. Work in this direction is in progress. In this paper we limit 
ourselves to give in Fig.^jthe squared Fourier transform of some quasihole functions, 

Mk) = J d 3 re lk - r Mr)- (45) 

In the valence states (lp for 16 and Id for 40 Ca) short range correlations slightly deplete |V>/i(fc)| 2 at large momenta 
with respect to the IPM. This behavior is in contrast with that of the total momentum distributions at large-fc, given 
in Fig.|| and showing a large enhancement due to the correlations. This discrepancy has been already observed f4S|| 
and it is confirmed by our approach. The effect of the correlations is more visible in the two low lying states given in 
the figure (Is for both nuclei). For instance, in the Is 16 case Jastrow correlations are effective at large momenta 
only, beyond the first IPM zero; instead, tensor correlations modify both the low and large momenta behaviors. The 
same effect is found in 40 Ca. 



V. SUMMARY AND CONCLUSIONS 



In this work we have calculated one-body density matrices, momentum distributions, natural orbits and quasi-hole 
states of 16 and 40 Ca using the FHNC/SOC resummation technique, which allows for using realistic interactions and 
state dependent correlations. The calculations have been done with the Argonne v' 8 two-nucleon potential plus the 
Urbana IX three-nucleon interaction, together with a set of single particle wave functions fixed to reproduce at best 
the empirical charge distributions of the two nuclei. The parameters of the correlation have been chosen to minimize 
the binding energies. Using these wave functions, we have investigated the role of the correlations on the quantities 
above mentioned. 

Our density and momentum distributions confirm some well known results. Short-range correlations have small 
effects on the density distributions and mainly around the center of the nucleus. On the contrary, the high momentum 
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tail of the momentum distribution is dominated by the correlations. We have pointed out that the tensor correlations 
enhance these tail by a factor 3-4 with respect to the results obtained with Jastrow correlations. 

The tensor part of the correlation is important in the calculation of the occupation probabilities of the natural 
orbits. The effect of reducing the occupation of the level below the Fermi surface and enhancing those which lies 
above is amplified by the tensor terms of the correlation. We found that the shape of the natural orbits below the 
Fermi surface is rather similar to that of the corresponding single particle wave functions. 

The natural orbits have also bcccn calculated within a lowest-order computational scheme. The agreement with 
the orbit below the Fermi surface is excellent. The lowest order calculation produces orbits above the Fermi surface 
with completely different shape with respect to those obtained by the full calculation. 

Tensor correlations play a relevant role also for the evaluation of the overlap functions and spectroscopic factors. The 
correlated overlap functions are close to the corresponding single particle wave functions if only Jastrow correlations 
are used. The inclusion of the tensor correlations strongly modifies their shapes. This behavior is also clear from 
the analysis of the spectroscopic factors. The depletion of a few percent, with respect to one, obtained with Jastrow 
correlations, becomes of about 10-15% for the valence levels and 30-45% for the deeply lying ones. 

In spite of this noticeable reduction, the FHNC/SOC approach in 16 is still unable to reproduce the empirical 
S P3/2 spectroscopic factor extracted from (e, e'p) reactions. A similar behavior was found in Ref. (5(| for nuclear 
matter, where the variational FHNC/SOC calculation of the one-hole strength, Z(e), around the Fermi level provided 
Z v (e ~ ej?) ~0.88, mostly due to tensor correlations. Second order perturbative corrections in a correlated basis, 
obtained by considering the contribution of two-hole one-particle, (2h — lp), correlated states, &2h-lp = G&2h-ip, 
were found to bring the strength to Zcbf{g ~ £f) ~ 0.70, so explaining almost half of the discrepancy with the 
empirical 208 Pb spectroscopic factor, Z( 208 Pb) ~ 0.5 — 0.6. The remaining part of the difference was attributed to the 
coupling of the single particle waves to the collective low- lying surface vibrations, not reproducible in infinite nuclear 
matter. We expect that in finite nuclei, the inclusion of correlated 2h — lp corrections can take into account also great 
part of the coupling with surface vibrations. 
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APPENDIX 



We present in this Appendix the FHNC/SOC expressions for the X n i and Af n i functions: 




(46) 



x {£f (r,rOQ"Vi)^r"Vi) ("P/PM^rr) + ^'(r.n)) 
+ J2 [fcrWi)^(r,ri)Cf (r x ) {-pf PM {r, n) + N^ c \r, r x )) 



p>2 




(47) 



+ £ [h^ir^h^ir^Ctfin) (-p? FM (r,r 1 ) + iV c c f i (r,r 1 )) 



p>2 
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<7^V,ri)QVi)^(r,ri) - JVg^r.n) - M^(r,n) A" A" 
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TABLE I. Occupation numbers of the nZ-th natural orbits for O and Ca in CBF, with the f§ and Jastrow correlation 
models. 



nl 


n ni (/ 6 ; 16 0) 


n ni (J; 16 0) 


n„i(/6; 4 °Ca) 


n ni (J; 40 C 


Is 


0.858 


0.960 


0.864 


0.952 


2s 


0.019 


0.005 


0.780 


0.962 


3s 


0.010 


0.002 


0.052 


0.002 


4s 


0.005 


0.001 


0.013 


0.001 


lp 


0.919 


0.980 


0.841 


0.949 


2p 


0.021 


0.004 


0.024 


0.009 


3p 


0.011 


0.003 


0.016 


0.006 


Id 


0.025 


0.006 


0.956 


0.983 


2d 


0.011 


0.003 


0.030 


0.007 


lid 


0.006 


0.001 


0.019 


0.006 
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TABLE II. The occupation numbers of the l-th natural orbits for le O and Ca in CBF, with three correlation models and 
in the Green function (GF) approach. 



corr. n s n p n d n/ 

*0 J 0.971 0.991 0.011 0.002 

U 0.977 0.988 0.015 0.004 

/e 0.899 0.966 0.047 0.006 

GF 0.936 0.951 0.022 0.007 

Ca J 1.932 0.975 1.003 0.022 

U 1-907 0.964 1.005 0.022 

fa 1.727 0.920 1.026 0.071 
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TABLE III. CBF spectroscopic factors for ie O and 40 Ca, with Jastrow (J) and spin-isospin correlations, with (fc) and 
without (/zi) tensor components. 

corr. Is lp Id 2s 

~ n7 J 098 098 

U 0.79 0.96 

fe 0.70 0.90 

40 Ca J 0.98 0.99 0.97 0.98 

U 0.71 0.76 0.96 0.97 

f 6 0.55 0.58 0.87 0.86 
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FIG. 1. Nuclear one body densities. The solid lines are the FHNC/SOC results with the fe correlation, the dot-dashed lines 
are the densities with the Jastrow correlation, the dashed lines are the IPM densities. 
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FIG. 2. FHNC/SOC momentum distributions in 1B and 40 Ca. Solid lines: fe model; dot-dashed: Jastrow; dashed: IPM. 
The squares are the VMC results of Ref. M . 
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FIG. 3. FHNC/SOC momentum distributions in ie O, 40 Ca and nuclear matter (NM). Solid line: NM f 6 model; dot-dashed: 
NM Jastrow model; thin solid line: 16 fa model; thin dashed line: 40 Ca fe model. 
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FIG. 7. Partial wave contributions to the one body densities in the natural orbits representation. The upper lines are the 
total densities. The remaining lines give v^^-pi(r,r). Solid lines: fe model; dashed: IPM. 
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FIG. 9. Squared transform of the quasihole wave function for O and Ca. Solid lines: fa model; dot-dashed: Jastrow; 
dashed: IPM. 
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